SIMULATING THE HOT X-RAY EMITTING 
GAS IN ELLIPTICAL GALAXIES 



Daisuke Kawata and Brad K. Gibson 

Centre for Astrophysics & Super computing, Swinburne University, Australia 
d kawata, bgibson@astro.swin.edu.au 

Abstract We study the chemo-dynamical evolution of elliptical galaxies and their 
hot X-ray emitting gas using high-resolution cosmological simulations. 
Our Tree N-body/SPH code includes a self-consistent treatment of ra- 
diative cooling, star formation, supernovae feedback, and chemical en- 
richment. We present a series of ACDM cosmological simulations which 
trace the spatial and temporal evolution of heavy element abundance 
patterns in both the stellar and gas components of galaxies. X-ray 
spectra of the hot gas are constructed via the use of the vmekal plasma 
model, and analysed using XSPEC with the XMM EPN response func- 
tion. Simulation end-products are quantitatively compared with the 
observational data in both the X-ray and optical regime. We find that 
radiative cooling is important to interpret the observed X-ray luminos- 
ity, temperature, and metallicity of the interstellar medium of elliptical 
galaxies. However, this cooled gas also leads to excessive star forma- 
tion at low redshift, and therefore results in underlying galactic stellar 
populations which are too blue with respect to observations. 

1. Introduction 

The hot X-ray emitting gas of elliptical galaxies represents an impor- 
tant interface between galaxies and the intergalactic medium (perhaps 
even the primary interface). The X-ray halos of ellipticals carry with 
them two fundamental mysteries: 

■ their X-ray luminosities are lower than that expected from an ex- 
trapolation of the cluster X-ray luminosity-temperature (Lx — Tk) 
relation (e.g. Matsushita et al. 2000). 

■ their X-ray metallicities are lower than that of the mean stellar 
iron abundance (the so-called "iron discrepancy" - e.g. Arimoto 
et al. 1997 - a "discrepancy" in the sense that the halo gas metal- 
licity was expected to exceed that of the stars, since it should bear 
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Figure 1 Dark matter density map of a portion of the 43 Mpc (comoving) 
simulation volume (upper panels), and predicted /-band image of the 
target galaxy (lower panels), over the redshift range z=3.0 to z=0.5. 

the pollution of the enrichment from earlier generation of stars - 
enrichment byproducts that were not locked up into subsequent 
stellar generations). 1 

Conversely, the optical properties of ellipticals appear less contentious! 
The Colour-Magnitude Relation (CMR) and Fundamental Plane pro- 
vide strong constraints for any elliptical galaxy formation paradigm. We 
present here our preliminary work aimed ultimately at the construction 
of successful self-consistent optical + X-ray cosmological chemodynam- 
ical simulations of elliptical galaxies. 

2. Methods 

In our simulations, the dynamics of collisionless dark matter and stars 
is calculated using a gravitational Tree N-body code, and the gas com- 
ponent is modeled using Smoothed Particle Hydrodynamics (SPH). We 
calculate radiative cooling, star formation, chemical enrichment, and su- 
pernovae (SNe) feedback, self-consistently, and take into account both 



1 X-ray iron abundances remain a controversial issue (c.f. Buotc & Fabian 1997), although the 
iron discrepancy appears to hold based upon recent high-resolution XMM RGS observations 
(Xu et al. 2002; Sakelliou et al. 2002). 
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Type la and Type II SNe. We assume that SNe feedback is released as 
thermal energy. Details of the code are presented in Kawata (2001) and 
Kawata & Gibson (2003, in prep). 

We have carried out a series of high-resolution simulations within the 
adopted standard ACDM cosmology (S7 =0.3, Ao=0.7, fJ b =0.019/i~ 2 , 
h=0.7, and 08=0.9). Gas dynamics and star formation are included 
only within the relevant high-resolution region (~12 Mpc at z=0); the 
surrounding low-resolution region (~43 Mpc) contributes to the high- 
resolution region only through gravity. The mass of individual gas par- 
ticles in the high-resolution region was 5.9 x 10 7 Mq. We next identi- 
fied an appropriate elliptical galaxy analog in the high-resolution region, 
which acts as the focus for this preliminary study. The total virial mass 
of this target galaxy is 2 x 10 13 M , similar in size to that of NGC 4472, 
a bright elliptical galaxy in the Virgo Cluster. The target galaxy is 
relatively isolated, with only a few low-mass satellites remaining at z=0. 

Figure 1 shows the morphological evolution of dark matter in the sim- 
ulation volume, and the evolution of the stellar component in a 200 kpc 
region centred on the target galaxy. The galaxy forms through con- 
ventional hierarchical clustering between redshifts z=3 and z=l; the 
morphology has not changed dramatically since z=\. Three different 
radiative cooling and SNe feedback models were considered: Model A is 
an adiabatic model (i.e. no cooling); Model B includes cooling and weak 
feedback; Model C mimics Model B, but incorporates stronger feedback 
(100 times larger thermal energy per supernova). 

For all the models, we examine both the resulting X-ray and optical 
properties, comparing them quantitatively with observation. The gas 
particles in our simulations carry with them knowledge of the density, 
temperature, and abundances of various heavy elements in their imme- 
diate vicinity. Using the XSPEC vmekal plasma model, we derive the 
X-ray spectra for each gas particle, and synthesize them within the as- 
sumed aperture (R~20 kpc). We next generate "fake" spectra with the 
response function of XMM EPN detector, assuming an exposure time 
(40 ks) and target galaxy distance (17 Mpc). Finally, our XPSEC fitting 
provides the X-ray weighted temperatures and abundances of various el- 
ements. Conversely, the simulated star particles each carry their own age 
and metallicity "tag" , which enables us to generate an optical-to-near in- 
frared spectral energy distribution for the target galaxy, when combined 
with our population synthesis code adopting simple stellar populations 
of Kodama & Arimoto (1997). 
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Figure 2 Left •panel: Comparison of the simulated and observed (crosses 
with error bars) Lx — Tx relations. The triangle/ circle/square indicates 
the predictions of Model A/B/C. The solid line indicates the extrapola- 
tion of the cluster Lx — Tx relation (Edge & Stewart 1991), while the 
dashed line represents the relation obtained from the adiabatic simula- 
tion of Muanwong et al. (2001). Right panel: Density vs radius (left) and 
density vs temperature (right) distributions of gas particles for Model A 
(upper) and Model B (lower). The sold curves separate the region where 
the cooling time is shorter (upper region) and longer than the Hubble 
time. 



3. Results and Conclusion 

The left panel of Figure 2 shows the predicted Lx — Tx relation for the 
three models at z=0; crosses with error bars represent the observational 
data from Matsushita et al. (2000). The adiabatic model (Model A) 
appears incompatible with the data due to its excessive luminosity and 
low temperature. The inclusion of radiative cooling leads to lower lu- 
minosities and higher temperatures - as a result, models with cooling 
(Models B and C) are (roughly) consistent with the Lx — Tx relation 
of the observed elliptical galaxies. These conclusions are consistent with 
the analysis of Muanwong et al. (2001). The right panel of Figure 2 
shows the effect of cooling more clearly. In the gas density versus tem- 
perature diagrams, the region above than the line corresponds to a pa- 
rameter space in which the cooling time is shorter than the Hubble time. 
Cooling ensures the gas within this region is cold (ie. non X-ray emit- 
ting), and of low density and high temperature, ensuring that ultimately 
radiative cooling drives the observed Lx — Tx relation. 

The left panel of Figure 3 compares the X-ray weighted iron abun- 
dance of our simulations with the observational data of Matsushita et al. 
(2000). As the adiabatic model (by construction) does not form any stars 
(not having any cooling!), we show the results only for Models B and 
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Figure 3 Left panel: Comparison of the simulated and observed (crosses 
with error bars) [Fe/H]x — Lx relations. The open circle/square shows 
the predictions of Model B/C; open stars show the predicted mean metal- 
licity of the stellar component for Models B and C. Right panel: Com- 
parison of the simulated CMRs (open circle/square for Model B/C) and 
that of the Coma cluster ellipticals (crosses). The dashed line shows 
the CMR fitted to the Coma Cluster galaxies. The solid circle/square 
demonstrates the colour and magnitude for Model B/C when the con- 
tribution from young stars (age <8 Gyr) is ignored. 

C. Both these models show lower gas-phase iron abundance, compared 
to their stellar abundance, consistent with the low iron abundances ob- 
served in the X-ray emitting gas of ellipticals. We find that a large 
fraction of iron ejected from stars is locked into future generation of 
stars. Stars preferentially enrich the gas in the central region, where 
cooling is efficient (right panel of Figure 2). The enriched gas can then 
cool easily and be incorporated into future generations of stars. Conse- 
quently, the hot gaseous halo has not been enriched efficiently, leading 
to a lower X-ray weighted iron abundance. 

In summary, our radiative cooling models explain the two X-ray "mys- 
teries" alluded to in Section 1. Having said that, any successful scenario 
must also explain the optical properties of the underlying stellar com- 
ponent. To this end, we examined the position of our simulated target 
galaxy in the observed Coma cluster CMR (Bower et al. 1992). We can 
see immediately that the colours of the resulting stellar components of 
both Models B and C are inconsistent with the data (being too blue). 
This inconsistency can be traced to an excessive population of young and 
intermediate age stars (<8 Gyr) which form from successively cooled gas, 
regardless of the strength of SNe feedback. If the contribution of these 
young stars was ignored, the resulting colours would match the observed 
CMR. Therefore, one exotic (if somewhat ad hoc) solution is to "hide" 
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these younger stars within a bottom-heavy initial mass function (IMF) 
such that they cannot be seen today even if they did exist (e.g. Fabian 
et al. 1982; Mathews & Brighenti 1999). Another (more plausible) pos- 
sibility is that extra heating sources, such as intermittent AGN activity, 
suppress star formation at low redshift. Before suggesting this is the true 
solution though, we must re-examine the predicted X-ray properties of 
the simulation end-products after introducing these additional heating 
sources; we will be pursuing this comparison in a future paper. 

Our cosmological chemodynamical code makes it possible to under- 
take quantitative comparisons between numerical simulations and obser- 
vational data in both the X-ray and optical regime with minimal assump- 
tions. We find that radiative cooling is required to explain the observed 
X-ray luminosity, temperature, and metallicity of elliptical galaxies. Un- 
fortunately, the resulting cooled gas also leads to unavoidable overpro- 
duction of young and intermediate age stellar populations, at odds with 
the observational constraints. Although a bottom-heavy IMF is one 
solution for this problem, extra heating by intermittent AGN activity 
seems to be more plausible (e.g. Brighenti & Mathews 2002); recent 
observations are consistent with this latter picture (e.g. Churazov et al. 
2001). 
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